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We give a direct path-integral calculation of the partition function for pure 3+1 dimensional 
U{N) Yang-Mills theory at large on a small S^, up to two-loop order in perturbation 
theory. From this, we calculate the one-loop shift in the Hagedorn/deconfinement tem- 
perature for the theory at small volume, finding that it increases (in units of the inverse 
sphere radius) as we go to larger coupling (larger volume) . Our results also allow us to read 
off the sum of one-loop anomalous dimensions for all operators with a given engineering 
dimension in planar Yang-Mills theory on IR^. As checks on our calculation, we repro- 
duce both the Hagedorn shift and some of the anomalous dimension sums by independent 



methods using the results of [hep-th/ 0412029| and |hep-th/ 0408178| . The success of our 



calculation provides a significant check of methods used in [hep-th/ 0502149] to establish a 
first order deconfinement transition for pure Yang- Mills theory on a small S^. 



1. Introduction 

Pure 3 + 1 dimensional Yang-Mills theory on has a discrete energy spectrum which, 
thanks to asymptotic freedom, can be computed perturbatively for small volume. This 
information is encoded in the thermal partition function, which can be used to evaluate 
thermodynamic functions and investigate the phase structure of the theory. In this note, 
we explicitly compute the partition function at two-loop order in perturbation theory for 
the U{N) Yang-Mills theory in the 't Hooft large limit [|l|, by a direct evaluation of the 
Euclidean path integral on x S^, where the radius of is Rss <^ 1/Aqcd and the 
circumference of the circle is the inverse temperature /3. 

Our results for the partition function may be expanded in powers of the dimensionless 
variable x = e~^^^s3 as 

oo 

Z2ioop{x) = ^x"(an + A6„ln(x) + 0(A2)), (1.1) 

n=2 

where a„ gives the number of states with energy n/Rsz in the free theory and 6^ A gives 
the sum of the order A perturbative energy corrections for all states with energy n/-Rg3 (A 
is the 't Hooft coupling Qym^)- Equivalently, bnX gives the sum of one-loop anomalous 
dimensions for all dimension n operators in the planar Euclidean Yang-Mills theory on M^. 
Using spin chain techniques, the matrix of one-loop anomalous dimensions has previously 
been computed in |^], so we can in principle derive the b^s from those results also. We 
find agreement between the two methods for all cases that we have checked explicitly. 

As discussed in [^^, large U{N) (or SU{N)) Yang-Mills theory in the limit of 
small volume has a Hagedorn density of states, manifested as a divergence of the partition 
function at a critical temperature given by Tc,o = [Rs^ ln(2-fy3)]-i at A = 0. At large (but 
not infinite) A^, this Hagedorn divergence signals a deconfinement phase transition across 
which the free energy jumps from order one to order A^^. As a physical application of our 
results, we can evaluate the leading perturbative correction to this critical temperatureB 
by determining the point at which our two-loop partition function diverges. We find that 
it increases (in units of the inverse radius) as the coupling increases i: 

T^Rss = T,,oi?53 ■ (l + + ^(^')) • (1-2) 

^ Note that beyond first order in perturbation theory, we expect the phase transition temper- 
ature to be below the Hagedorn temperature Q, but to order A, they are the same. 

^ Since the coupHng runs, we should clarify that the relevant coupling is the one at the scale 
l/Rgz, so an increase in the coupling is the same as an increase in RgaAqcD- 
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This is consistent with the conjecture that T^Rs^ is monotonic in the radius, since at large 
radius we expect Tc to approach a constant of order Aqcd^ so T^Rs^ is an increasing 
function of i^^a. As another check of our results, we provide an independent calculation of 
the shift in using a formula in for the shift in the Hagedorn temperature in terms of 
the matrix of anomalous dimensions, which we can get from |^]. Again, we find that the 
two methods give the same result. 

While our results are interesting in their own right, one of the main motivations for 
this work was to provide a check of the calculation and regularization methods used in 
1]^ to establish that the planar Yang-Mills theory on a small undergoes a first-order 
deconfinement transition as the temperature is increased. That calculation involved a 
number of novel features related to the spherical space: sums over angular momenta instead 
of spatial momentum integrals, integrals over three- vector or scalar spherical harmonics on 
appearing in the vertices, and a novel non-gauge-invariant cutoff regularization scheme 
used together with counterterms chosen to restore gauge invariance in physical calculations. 
While the calculation in ||^ had internal consistency checks, such as the cancellation of all 
logarithmic divergences among the diagrams, we did not have any independent method of 
verifying the final numerical result (the sign of which determined the order of the phase 
transition) . Since the calculational setup and many of the steps are identical in the present 
calculation,! we view the successful matching of our results here with those from other 
methods (involving only standard calculations on IR^) as a satisfying check of the validity 
of our formalismB. 

The structure of the paper is as follows: in section 2, we provide the basic setup for 
our calculation and outline the two-loop diagrams that will be needed in order to com- 
pute the partition function. These are divergent, so we need in addition certain one- loop 
counterterm diagrams. In section 3, we review our regularization scheme and calculate 
the necessary counterterms, including a new curvature-dependent counterterm not present 
for fiat-space calculations. In section 4, we evaluate the regularized two-loop diagrams 
together with the one-loop counterterm diagrams in order to get our final result for the 



^ The calculation in [|| was more specific in that the necessary information was contained in 
specific two and three-loop contributions to the effective action for the constant mode of Aq on 



X after the path integral was performed over all other (massive) degrees of freedom 
^ In fact, the computation described in this paper helped us to find a small mistake 
original version of the calculation of M]. This mistake did not affect the main result of M, 



partition function and, subsequently, for the one-loop correction to the Hagedorn temper- 
ature. We verify that our results are independent of the regulator used. In section 5, we 
verify our results by independent calculations of the Hagedorn shift and of the first several 
terms in the expansion (|1.1| ) using results from and [||. 

2. Setup 

The basic setup for our calculation is identical to that in [H, so we give only a brief 
summary here, referring the reader to section 2 of that paper and section 4 of for more 
details. 

We would like to calculate the thermal partition function 

Z = Y, e"^""^ = e-^^^'"'^' = (2-1) 

it i 

for pure U{N) Yang-Mills theory on with radius -R53, where /3 is the inverse temper- 
ature, Ai are the dimensions of the local operators in the theory, and Ei are the energy 
levels. This is given by the Euclidean path integral on xS^, with the circle (compactified 
Euclidean time) chosen to have circumference P, of the Euclidean Yang-Mills action 

Seuc = \jytj d^'x tT{F^^F^^) . (2.2) 

We use a normalization in which the Yang-Mills coupling Qym appears in the interaction 
terms in the action. To make the path integral well-defined, we need to fix a gauge 
and include the appropriate Fadeev- Popov determinants. For our calculation on S'^, it is 
convenient to use the gauge 

^^A' = 0, dt [ Ao = 0, (2.3) 

where i runs over the spatial directions. The latter condition implies that the zero-mode 
of Aq on is constant in time. We denote this mode by a. 

Apart from a, which has no quadratic term in the Yang-Mills action ( |2.2| ), all other 
modes are massive. As such, it is convenient to first integrate them out to obtain an 
effective action for a. As argued in [^, this effective action can only depend on the unitary 
matrix U = e*^" (the Wilson line of the gauge field around the thermal circle, averaged 
over the sphere). Further, the effect of the Fadeev- Popov determinant associated with the 
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second gauge-fixing condition in (^]^) is to convert the measure [da] in the path integral 
to the Haar measure [dU]. Thus, we obtain 

[(/;7]e-Se//(f/) (2.4) 

where 



g-s.//(c/) _ j [c;A'][(ic][(ic]e-^^"=("^''")-^^^(^'"). (2.5) 

Here, [dA'] denotes the measure for the gauge fields excluding the zero mode of Aq. The 
fields c and c are ghosts and 

Sfp = - f dt f tx{dic&c) (2.6) 

is the ghost action (where is a covariant derivative) , introduced to represent the Fadeev- 
Popov determinant associated with the first gauge-fixing condition in (2.3). 



2.1. Vertices and propagators 

To perform the calculation, we expand all fields into modes on S"^ (we set the radius 
of 5"^ to be i?53 =1 from here on as it can always be reinstated by dimensional analysis), 

Ao(t,^) = ^a°(t)^"(^); 

a 

A,{t,6) = Y,APmhO): (2.7) 

a 

c(t,^) = 5^c"(t)5"(^), 

a 

where 5"" and Vj^ are scalar and vector spherical harmonics on 5""^, defined in appendix 
B {a and /3 represent angular momentum quantum numbers). This leads to the action 
in appendix A. Because the fields a" and c" appear only quadratically, we can integrate 
them out first to get an effective action involving only A^ and a. The final result includes 
a quadratic action 

S2= I dt tr (^A^i-Dl + + 1)')^") , (2.8) 



a cubic interaction 



^3 = QYM j rfttr(zA"A^A^e«(i„ + 1)^;"^^), (2.9) 



and quartic interactions 



where we have defined D-^ = dt — i[a, ■], and D'^^t and E"-^"* are integrals over products 
of three spherical harmonics that are defined in appendix B. Sums over all indices are 
implied. 

There are in addition higher order vertices arising (like the second line in ( |2.10|) ) from 
integrating out a and c, but these do not enter in our two-loop calculation. We have 
also suppressed a set of vertices proportional to 5(0) which serve to precisely cancel terms 
proportional to 5(0) that arise in contracting two I^t-A's at equal times. As explained in 
1^, we can simply work with the set of interactions above, ignoring any 5(0) terms that 
arise.i 

The propagator for A" that follows from ( |2.8| ) is 



{K,miait')) = 5-~^^ff\t - t\ a), (2.11) 
where A is a periodic function of time given for < t < /3 by 

^ 2{j + 1) \l-e'»f}e-{j+i)f3 - I _ eza/3eO-+i)/3 J ' (^.12) 

Here, a is shorthand for a ® 1 — 1 ® a, and a term a"^ ® in the expansion of A should 
be understood to carry indices (a'")"'^(a"^)^^ in ( pTl| )- 

2.2. Regularization and counterterms 

The two-loop diagrams that we are required to evaluate contain divergences, as we 
should expect from gauge theory in four dimensions. To deal with these, we could apply 
dimensional regularization, evaluating the partition function on x x IR"^, and using 
a modified minimal subtraction procedure to obtain finite and gauge-independent results. 



^ These 5(0) 's are related to our choice of Coulomb gauge, which leads to singular propagators 
for Aq and the ghosts. In appendix D, we explain a more well-defined way of performing Coulomb- 
gauge calculations which avoids any 5(0) terms in the calculation but nevertheless gives the same 
results, justifying our naive cancellations. 
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In practice, to overcome technical difficulties with this regularization scheme, we will use 
the regularization scheme of involving momentum cutoffs, which requires adding to the 
action counterterms designed to yield results that are completely equivalent to those of 
dimensional regularization.i This is reviewed in detail in section 3. For our calculation, 
the counterterm vertices that contribute arei 

= xl'dt tr (Af ^^^^ {Zo - Z,d] - Z^D^Ar^"^^) 

- A / dttr {A"{Zo + Ziij^ + If - Z2Dl)A^) (2.13) 
A /"^ 

where A = Qym-^ regulator dependent constants that we will determine in 

section 3. 

2.3. One loop result 

The evaluation of the partition function at one-loop order has been carried out in 0] . 
The result is 



where 



and we define 



^i-ioop = J [rfC/]e-^?-'°°p(^), (2.14) 

oo 

Sf-ioo^U) = -Y, -^(x-)tr(t/-)tr(t/tn), (2.15) 



n 

n=l 



- 2x^ 

which is the single mode partition function for a free vector field on S^. The unitary matrix 
integral can be evaluated explicitly at large to give @] 



-^i-ioop — e 



-/3_Fi„ioop 

1 — z{x 



oo ^ 

n 1 — r^- (2.17) 



More precisely, the results are equivalent to dimensional regularization in a gauge where 3 + d 



components of the gauge field participate in the Coulomb gauge condition (2^). This is related 
to the more general approach of split-dimensional regularization Pj8|,p|,p!0[ , in which the number 
of degrees of freedom that do not participate in the Coulomb gauge condition is also varied. 

As indicated, it is only the SU{N) part of the gauge field Af^^'^^ = Ai — j^tv{Ai) that can 
be present in the counterterms since the U{1) part of the theory is free. 
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The function z increases monotonically from to oo as x increases from to 1 (i.e. as the 
temperature increases from to oo), so the n = 1 term in the product leads to a divergence 
as the temperature is increased to the critical value Xc,o such that 2(xc,o) = 1, or 

7^0,0^53 = (ln(2 + VS))"^ ^ 0.75933. (2.18) 

This is the Hagedorn temperature of the large free theory (and of the interacting theory 
in the small volume limit). For higher temperatures (at finite A^) there is a different saddle 
point dominating the path integral, and Z behaves as exp {—(3N'^f(T)). 

2.4- Two loop calculation 

At two-loop order, the partition function is given by 

^2-ioops = J [dt/]e-^--p(^)-5f-!>oo.(^), (2.19) 



where 

e = (e ^'"*)2-ioop 

/ 1 ^ (2.20) 



2 

Here, the expectation values are evaluated in the free theory with fixed a. 

The correlators in ( |2.20| ) contribute to the partition function in two different ways 



TlJ|. First, as shown in H, the planar diagrams give a contribution of the form 



S;f = PglMNY.fnixMU'^MU^n+Pgl^ ^ /^^(x){tr([/-)tr(t/"^)tr(C/-— )+c.c.} 

n=l n>m>0 

(2.21) 

The three-trace terms do not contribute to the partition function at order A, so they can be 
ignored. On the other hand, the double-trace terms modify the Gaussian integral ( |2.14| ) 
(the path integral is Gaussian in the variables Un = tr{U^)/N), and result in order A 
corrections to the denominators in ( |2.17| ). 

Since (|2.17|) is actually the sub-leading contribution to the large free energy (the 
leading 0{N'^) contribution, coming from the action evaluated on the saddle-point, van- 
ishes here), there are also contributions at the same order arising from non-planar two- 
loop diagrams. These are independent of U (they have a single index loop) and give a 
temperature-dependent prefactor to the infinite product in ( |2.17| ), 

Sf,^/39'YMNF-^{x). (2.22) 
7 



Thus, in terms of the functions /n(x) and defined in (|2.21|) and (|2.22| ), the 

final result for the two-loop partition function is 



2; = e-Wr(^) y" \dV\^^ ^(2(x")-A/3n/„(x))tr([/")tr(t/t") + . 



00 



(2.23) 



n=l 



Expressed in terms of the correction to the free-energy, we have 

5F = A(^F-(x)-ff:^^^j. (2.24) 

From (|2.23| ), we see that the corrected partition function will diverge when 

l-2(x) + A/3/i(a;) =0, (2.25) 
so we find that the critical value of x shifts by 

5xe = A/3,,o4^4 (2.26) 

or, equivalently, the critical temperature shifts by 

= AT,o-M^ . (2.27) 

It remains to evaluate /n(a^) and F^^ix) by evaluating the planar and non-planar two- loop 
diagrams plus one-loop counterterm diagrams. We do this in section 4, but first we must 
discuss our regularization procedure and determine the necessary counterterms. 



3. Regularization and counterterms 

The regularization procedure that we use was described in detail in ^j, so we only 
summarize it briefiy here. The central idea is to cut off angular momentum sums in 
such a way that our calculations are as simple as possible. Such a scheme will in general 
break gauge and Lorentz invariance but, by choosing the right set of local counterterms 
in the action at the cutoff scale M, we can ensure that both are restored in the theory 
far below the cutoff. Specifically, we will choose counterterms so that the results for any 
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correlator match precisely, in the limit where the cutoff is removed, with those obtained 
using dimensional regularization. 

Though it would be conceptually simpler to perform our calculations directly in di- 
mensional regularization, this presents technical challenges for the required diagrams that 
we have not been able to overcome. On the other hand, determining the appropriate set of 
counterterms to ensure the equivalence of a more general cutoff scheme with dimensional 
regularization requires, at least to low orders in perturbation theory, the evaluation of only 
a few simple diagrams using both methods and a comparison of the results. 

The regulator that we employ requires the insertion of a damping factor R{ a/p^/M^) 
for each Ai propagator, where p'^ = piPi {i = 1,2,3) is the magnitude of the spatial 
momentum and R is a function satisfying -R(O) = 1, -R'(O) = 0, and R{x — > oo) = i. The 
regulator preserves rotational invariance, but breaks Lorentz and gauge-invariance, so we 
should include all possible counterterms which are local, rotationally invariant, and have 
dimension four or less. Since all counterterms will be at least of order A, only quadratic 
counterterms (giving rise to one-loop diagrams) can contribute to the partition function 
at order A. Moreover, only the Ai propagators are temperature-dependent with our choice 
of gauge. Thus, only the counterterms appearing in ( p.l3| ) can contribute to our resulli. 

In order to determine the constants Zq, Z\, and Z2, we need to evaluate at least 
three simple one- loop correlators to which Zq, Z\ and Z^ contribute in different linear 
combinations, and choose the Z's so that the results match with the same correlators 
evaluated in dimensional regularization, with a minimal subtraction scheme to be described 
below. 

It is important to note that for our calculation on the sphere, the Zq counterterm 
actually combines three separate local covariant structures, 

y" d^'x^ (Z(f'"*tr(A,AO H- 4mr(A,A,) + Z^'7^,,tr(A,A,)) , (3.1) 

the last two involving the Ricci scalar and Ricci tensor built from the metric. If Zg or Zq 
are nonzero, the \x{^AiAi) counterterm on the sphere will be different than the one in flat 

^ Note that we do not use a regulator for Aq or ghost lines, which we have already integrated 
out explicitly. As described in appendix D, the cancelling 5(0) terms that appear in integrating 
out these fields are not physical UV divergences and may be regulated by choosing a more general 
gauge. 

^ Note that a possible term of the form diAidjAj will not contribute due to the Coulomb gauge 
choice. 
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space, with coefficients proportional to 1/ R^a in addition to tlie flat space coefficient Zg 
proportional to M^. Thus, it is important that the calculation used to determine Zq be 
performed on the sphere, as we will do in section 3.2. 

On the other hand, the counterterms involving Zi and Z2 are already dimension 4 
operators, so there is no allowed structure involving the spatial curvature that reduces to 
these. Consequently, it is enough to study fiat-space correlators to determine Zi and Z2, 
and we turn to this presently. 

3.1. Curvature-independent counterterms 

The Zi and Z2 counterterms were already determined in [|5|, but for completeness, 
we review the essential parts of that calculation here. As we suggested above, it is sim- 
plest to determine the curvature-independent counterterms by a flat space calculation. To 
determine Zi, we will compute 

(A,(0,p)A,(0,-p))jr4 (3.2) 

i.e. the term proportional to p'^dij in the IPI two-point function of Aj, while to determine 
^2, we will compute 

(A,(a;,0)A,(-cc;,0))i?^ (3.3) 

with CO that component of the momentum which does not participate in the Coulomb 
gauge constraint. In each case, we compute the result (to order A) using our regulator and 
the counterterm with undetermined coefficient, then repeat the calculation in dimensional 
regularization with a modified minimal subtraction scheme. We finally determine the 
counterterm by demanding that the two calculations agree. 

Dimensional regularization 

We will begin with the calculation in dimensional regularization. We generalize our 
Coulomb gauge by assuming that 3 + d components of the gauge field participate in the 
Coulomb gauge condition. With this choice, we may write the quadratic action as 

^ ^ laM + laM ^ m,c] . (3.4) 

The interaction terms for Euclidean Yang-Mills theory on IR^"*"^ in the Coulomb gauge 
include a cubic action 

S3 = 9YM j d'^^'^xtr^^-id,Aj[A,,Aj]-iA,[AQ,A,]+idiAo[Ao,Ai]-idic[Ai,c]^ (3.5) 
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and quartic terms 



(3.6) 



From the quadratic action, we may derive propagators (suppressing the color indices) 



A;2(z/2 + A;2) 



(Ao(z/, /c)Ao(-z/,-/c)) 
(c(z/, k)c{-iy, -k)) 



1 

k^' 
1 



(3.7) 





SEla 



SElb 





SEld 



SEle 



SElc 



SElf 



Figure 1: Diagrams that contribute to the two-point function (AiAj). 
Dashed (arrowed) hues denote Aq (ghost) propagators, and sohd hues denote 
Ai propagators. 



The quantities ( |3.2D and ( |3.3| ) are both obtained from the IPI two-point function of 
the gauge field. At one- loop, this gets contributions from the six diagrams of figure 1. 
Note, however, that diagrams SEle and SElf, which involve only quartic vertices, do not 
depend on the external momenta and hence do not contribute to (^3) or (|3.3|) . In addition, 
the Aq loop diagram SElc cancels with the ghost loop diagram SEld0 so we need only 
focus our attention on the first two, SEla and SElb. For these, we find 



Actually, these two diagrams only cancel up to a 5(0) divergence of the usual sort associated 
to ^0 lines. As discussed in Appendix D, this divergence is completely unphysical and is not 
important for any of our results. 
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-]^{A,{u^,p)A^{-u,-p))Ye{^ 

-j^-^ {-A,, [(p + fc) ■ A ■ b + k)\ + [b + fc) ■ A ■ A],(2/c - p), 

+ /c) ■ A ■ A],(2A; -p), + [(2p - A;) ■ A],[(p + /c) ■ A], + A^AkMp - 2k) , 

(3.8) 

where A = A{iy, k) and A = A(a; — u,p — k), and 

--^{M-^P)M--^-p))l'A = j (3.9) 

It is straightforward to insert exphcit expressions for the propagators and expand to order 
p2 or u'^. In the latter case, we findllll 

~{A,{u;, 0)A,(-a;, 0))^?^ = ^'S,,^ {-/o,i(0) + 2/0,3(0) - 8/1,4(0)} , (3.10) 

where we define 

d^+'^k dv zy2m^2n-2m-2 



{27t)^+<^ {k^ + a^){u^ + k^)'' 

r(n-m-^)r(m + |)r(i + f)r(-f; 

5 + d 



24+V(7r)^ r(n)r(f + f) 

1 r(n-m-^)r(m + ^) / //xx 
47r3 Tin) \ \aJ 



(3.11) 
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Here, /U is the regularization scale introduced, as usual, to keep /m,n(a) dimensionless as d is 
varied. In the last line, we have set d = —e and expanded for small e. The integrals for a = 
in (|3.10| ) contain infrared divergences in addition to the UV divergences, so rather than 
comparing these directly between the two schemes, we write /(O) = /(a) + {/(O) — /(a)}, 
and, noting that the expression in curly brackets contains no UV divergence and must 
agree between the two schemes, focus on the /(a) term and compare this between the two 
schemes to determine the counterterms. 

For our calculations, it is convenient to choose a minimal subtraction scheme which 
sets the combination in square brackets in ( |3.11|) to zero. With this choice, the only non- 
vanishing contributions to (|3.10| ) are terms proportional to ln(|U/a) and terms for which 



We have used the fact that kikj inside the integral may be replaced by k^5ij/{3 + d) if the 
external momentum p is set to zero. 
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an e in the expansion of the (i- dependent coefficients multiphes the ^ part of some Im,n- 
The result for fE73|) is thus 



— {Mu, 0)M-u, 0))i?^ = - ^ In (^) + {UV finite}) , (3.12) 

where the UV finite term denotes the difference between the original integral and the IR 
regulated versionlll. Following the same steps, we find that 

(3.13) 

Cutoffs and Counterterms 

We now reevaluate the correlators using the regulated momentum integrals together 
with counterterms, choosing the counterterm coefficients so that the results match with 
those of ( ^.12| ) and ( ^.13| ). Our regularization scheme employs a cutoff only for the Ai 
lines, so the regulated expressions for diagrams SEla and SElb follow by setting d = and 



replacing Aij(z/, /c) /S.ij{y,k)R{k/M) in (|3.8| ) and (|3.9|). Our results may be expressed 
in terms of the following basic integrals: 



Cn^^l dqqR^iq), (3-14) 







dqqR{q)R'\q). 







Again, we start with ( p.3|) , for which we can set p = 0. In this case, the contributions 
from ( p.8| ) and ( p.9|) now include regulator factors of R^{k/M) and R{k/M), respectively. 
Starting from (|3.10| ) with (i = 0, we can evaluate the integrals over u to get 



1 r(n-m-i)r(m+l) f d^k 1 
"^'"^ ^ 2TV r(n) J (27r)3 fc(fc2 + a2) ^^"'^^ 



^2 Note that this term will contain a dependence on the IR cutoff ln(a), which exactly cancels 
that explicitly written. 
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and then insert the regulators. Including also the counterterm contribution, we find 



-'^d^^-^ln^ +^ln 



(3.16) 

For the other correlator (|3.2| ), we again need a regulator factor R{k/M) for the contribution 
{ ^ ) from SElb, and this time the factor R{k/M)R{\p - k\/M) for the contribution (U) 
from SEla. It is important to take the expansion of R{\p — k\/M) in powers of p into 
account in order to correctly obtain the terms proportional to p'^Sij in ( |3.8| ). The effect of 
this is trivial for logarithmic divergences as only the first term, R{k/M), contributes there. 
However, the expression ( p.8|) has a quadratic divergence, for which subleading terms in 
the expansion become important. Working everything out carefully, we find 

1 ( r rPk ^ /^ 2 

--{A{0.p)MO^-p)y/l=p%,}^Z^ + J (2,)3fc(fc2 + ^2) [-^Rik/M)--Rik/Mr 

1 , fAiM\ 9 , fA2M\ F2 1 



-4 



p^d,,{Zi + ^^\n( '-^ ) -:^ln 



IOtt^ \ a J 407r^ \ a J 15 40n 

(3.17) 

where the first line contains the contribution from as well as the purely logarithmically 
divergent terms in (3^) and the second line contains the only term in ( p.8|) for which 
R{\p — k\/M) must be expanded in p beyond the leading order. Finally, demanding that 
( p7| ) and ( pJ6|) agree with ( p^ and (|3J2D , we must have 



1 

87r2 


In 1 


K ^ J 




—F2 + 
15 


1 


In 1 






1 


87r2 


K ^ J 


487r2 



487r2 407r2 \Al 
1 , M2 



(3.18) 



3.2. Curvature dependent counterterm 

As we discussed above, the tr{AiAi) counterterm has a coefficient which depends on 
the spatial curvature, so we need to determine Zq by a direct calculation on 5"^, though 
we can work on x IR^ rather than x . Thus, to find Zq we compute the two-point 
function of a mode of Ai on the sphere, again demanding that a calculation using regulator 
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functions and counterterms matches with the result using dimensional regularization (this 
time gauge theory on IR^"'"'^ x 5"^). Specifically, we will compute 

R= J2 (^'"^'"'"'''(^ = 0)A^=^'"''"''''(c^ = 0))ipj , (3.19) 

the IPI two-point function of the lowest total angular momentum mode of Ai, summed 
over polarizations e and angular momentum states m,m'. We begin with the calculation 
in dimensional regularization. 

Dimensional regularization 

To define the dimensionally regularized theory, we choose 

Seuc = J d''+^xcPytT {^^FuFij^ , (3.20) 

where y are the coordinates on 5"^ and x are the coordinates on IR'^"'"^ with d = — e, and 
for simplicity we take the regularization scale to be unity (in units of Rga), 

H = l. (3.21) 

We denote by Aa the components of the gauge field in the d directions, while Aq and Ai 
denote as before the components of the gauge field in the time and sphere directions. We 
work in the Coulomb gauge, now extended to involve d + 3 components^! 

d,A, + daAa = . (3.22) 

Since the divergence of Ai is no longer set to zero, the expansion of modes on the sphere 
now becomes 

a 
a 



For this part of the calculation, it is actually much simpler to choose a gauge for which only 
the sphere components of the gauge field participate in the Coulomb gauge condition, but we 
must make the present choice to be consistent with the conventions of Q and of the previous 
subsection, used to compute Zi and Z2. 
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Note that the Coulomb gauge condition determines the coefficients of VS in the expansion 
of Ai in terms of the Aa modes. With these expansions, we find that the quadratic action 
may be written as 



+ la^-dl+Mja + 2))a'^ (3.24) 

From these, we determine the propagators to be 
(a'^(a;, k)a^{u, k)) = {27r/+''S{uj + u)5'^{k + k) ^'^'^ 



{A'^iuj, k)Af^{u, k)) = {27rf +'^5{uj + u)5'^{k + k) 



c^2 + A;2 + (j„ + l)2' 



(3.25) 

There are now additional interaction terms involving A. For our calculation of the two- 
point function of Ai, we will need the cubic terms with at least one power of this field, and 
the quartic terms with either two or four powers of A. In addition to those terms which 
appear for d = 0, listed in (A. 2) and (A. 3), the new terms of this form are: 

C/3a7 



CAaA=9YMtA-idodaAZ[a^,A'^] . - idoA^ia^ , d^Al] . .. 

J^AAA =9YMtr iiA^A^AiC^^^ - i[d,At, daA^]Al^ 



-iK, dad,A^,]Al - iA-daA%Al ^l^J' C^--^ , 

Jp (j^ + 2) 313 [313 + 2) (j^ + 2) J 

(3.26) 

and 

>C1 =^yMtr [-\[A'',,Af'][Al,A']DP'^B'^^~>^ - ^-[A^,daAfi][A\d^Ai]D'^'^^m^ 
--[A-,daA^,][d,AlA']C>^-^C>^'^ - -[A-,A^][daAZ,d,At]-^ 



2' ^—an-'^-b^ J- - 2^ ' ^^-'^•---''•-''^^•^(^•^ + 2)^,(^5 + 2);' 

(3.27) 
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SElg 



SElh SEli SElj 

Figure 2: New diagrams with Aa propagators that contribute to the two- 
point function (AiAj). Dotted (dashed) hues denote Aa (a) propagators. 

where S, C and B are defined in appendix B. 

We are now ready to calculate ( |3.iy| ), setting the external momentum in the d di- 
rections to zero also. In addition to the diagrams of figure 1, the new interaction terms 
( p.26|) and (|3.27|) give rise to the diagrams of figure 2. Note that because the latter con- 



tain internal A lines, they yield contributions to the two point function only through the 
multiplication of 1/e poles (from sums which are logarithmically divergent at e = 0) and 
prefactors proportional to e. 

For the n'th diagram, we define 11"^ to be the correlator for two arbitrary modes 
A" and A^. This has a contribution proportional to 5(0), which we denote by 5n (these 
cancel between the diagrams, and the naive cancellation may be justified by the method 
in appendix D), and a finite contribution. After setting (3 = a, choosing = 1 and 
summing over the external m's and e's, we denote the finite contribution by Gn- Some 
useful formulae may be found in appendix C. 

We begin with the diagrams involving quartic vertices. Diagram SEle gives 

n... ^ 2(I,--I,«- - D"-!,-^) / ff.,.^,.;^,^^^. . (3.28) 
From this, we find 



(3.29) 
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We turn next to diagram SElf 

This gives (at order e°) a 6(0) term 

SsEif = -25(0)(I?°^^I?^^^ + ■ ,.\ r7^aA^^/3A) / (3.31) 
and no remaining finite contribution, 

GsEif = 0. (3.32) 

We move now to diagram SElg, which receives contributions from the four quartic 
interaction terms in ( |3.27| ). The first such term contributes 

n...,, ^ -2Z,-i,-^ / |^_,_^,J_^„ - ,,^^). (3.33, 

We find that the two terms in the integral actuaUy cancel each other, so the net result is 

GsEig,i = . (3.34) 
For the second term in ( |3.27|) we find 

Use.,,.- -2D B J (2.)^+i + p +,-^(,-^ + 2))(P +,,(,, + 2)) ' ^''''^ 

To evaluate this, we perform the integral over lo followed by the k integral. The resulting 
summand may be expanded in powers of and since we find an overall factor of d, the 
only contributing term is 

-ii:^-T^-|.C(l + .)-i+OW. (3.36) 
Thus, this diagram gives a net contribution 

3 

GsElg,2 = (3.37) 

From the third term in ( |3.27| ), we find 

''''''' ~ + J i27TY+Hu^' + k'+h{h + mk'+j,{jj + 2))- ^ 
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The evaluation of this is similar to the previous diagram, and we find a net contribution 
of 

GsEig,3 = (3-39) 



Finally, the last term in ( 3.27 ) gives no contribution. 

^5^13,4 = (3.40) 
We now move to the cubic diagrams. Diagram SEla gives 
n _ 0^075^/375 { dud'^k 1 1 ('iA^^ 

In this case, choosing = 1 forces js = jj by the triangle inequality. We find 



GsEia = ^-^nl - ^)(C(-1 -d)+ 4C(1 -d)- 5C(3 - d)) 

8,1 1 , , ^ 1 ^ 47 10^,^, , 
TT^ e 2 2 DTT^ TT^ 



(3.42) 



Diagram SElb gives 



dud'^k ( 1 /c2 + (jp + l)^ 



(27r)^+i + JtO't + 2) (^2 + J70'7 + 2))(a;2 + p + (_^-^ + 1)2^ 

(3.43) 

The first term here gives a 5(0) term 

5sEib = 2S{0)D^P^Df'^^--^—— (3.44) 

J7U7 + 2) 

which cancels the first term in SsEif, while the second term gives a finite contribution 



GsEU = '-^{- + 1 M^r) + J7) + A + 0{e). (3.45) 

TT^ e 2 2 TT^ 



For diagram SElc, we find 



This gives only a 5(0) contribution, 

Ss... ^ -,mC^'C'^>-^-^±^^ (3.47) 
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with no remainder, 

GsEic = 0. (3.48) 
The ghost loop diagram SEld also gives only a 5{0) term 

Ss.u = 2.(0)C-"C'^^-y-^^i^y-^ (3.49) 

which, together with Sseic, cancels the second term in SsEif- 

The remaining diagrams contain internal A lines, so for each of these, it is only- 
necessary to isolate the logarithmically divergent term in the sum. From diagram SEli, 
we find 

_ C^^^^C^ f dud'^k fc^ 

" " jx{jx + 2) J (27r)'^+i (F + + 2))(a;2 + P + + 2)) ' ^ 

where we have subtracted a term which gives a 5{0) from the oj integral because it is 
multiplied by a A; integral which vanishes for e — > 0. This gives a finite term 

GsEii = -^^^ (3.51) 
For SElh, we have the contribution 



/ f dud k 

3x{3x + 2) J j^-^^ 



jx{jx + 2) jxUx+2) 



(^2 + + (^-^ + 1)2)(^2 + ^2 + j^^j^ + 2)) fc2 + j^^j^ + 2) • 

(3.52) 

This gives a finite contribution 

4 

GsEih = ■ (3.53) 

Finally, we find that diagram SElj gives a net contribution of 

dud'^k 1 



^SElj = -CP'^'C^'^P j 



i2nY+i (^2 + k^^3.{3. + 2))(c^2 + p + 2)) 

• [M + 2/c2(S, + - 2(Ep + E,)) 

+k\{B^ - 2E,){B^ - 2Ep) - 2BA - 2B,3Ep) 

-k\iBp - 2Ep)BaE^ + (Ba - 2E^)BpEp) + k^B^BpEpE^] , 

(3.54) 
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where we have defined 

p 1 „ _ j^j^ + 2) + jpjjp + 2) + {j^ + 1)^ 

^ " + JpOp + 2) ' - + 2)j,ij, + 2) • ^^-"'^ 

In evaluating this, things simphfy, since we take = jp = 1 and this forces jp = We 
have, as an intermediate step, 

^ ^ a(a + l)^(a + 2) f d'^k 1 / (l+jSfc>(a + 2) \ 

^""^^■-4- y (^(fc2 + a(a + 2))i\^^~'^+l F + a(a + 2) J 

(3.56) 



where 



a(a + 2) + 2 , , 



After computing the integrals, we find all terms are proportional to d. Expanding in powers 
of 1/a and keeping only l/a^~'^ terms in the summand, we find a net contribution of 

GsEij = --^- (3-58) 
Combining all terms, our total result is 



r -A 

^ dimreq o 



11 1 , . N 1' 



+ -^7+ A-^C(3). (3.59) 



Using the modified minimal subtraction scheme defined in §3.1, the quantity in square 
brackets is set to zero, so we find finally that 

6 3 10 



Gdimreg = "7 + TT^ jC(3). (3.60) 



Cutoffs and Counterterms 

We now repeat the calculation using a regulator function -R(a/— V^/M^). Our scheme 
applies this only to the Ai lines, for which — gives [j + 1)^ for the mode with total 
angular momentum quantum number j . The expressions Hi for the diagrams are the same 
as above, so we simply set d = and insert the regulator functions. We express our 
results in terms of the basic integrals defined in (|3.14|) , using the Euler-McLaurin formula 
to compare infinite sums with integrals. 

For diagram SEle, we find 

a=l a=l ^ ' 

(3.61) 
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(3.63) 



For diagram SEla, we find 

G'sE^a = 4 E -R'O + 4 E -R'O - ^ E 4 

a=l a=l a=l (3_g2) 

Diagram SElb contains a 5(0) piece that is cancelled by the remaining diagrams. This 
leaves the result 

a=3 a=2 

3 2 a(a^ ^ ■*■) d/- ^ \ 
= - - ^ (a2 -4) 

23 2 ■^S^ , a 6 -'■ D,' \ 

TT^ \ fi J 37r^ TT^ 

Combining all terms so far, we have 

\ fj. / TC^ \ ^ J DTT^ TT^ TT^ 

(3.64) 

X 

Figure 3: Counterterm contribution to (AiAj). 

We have in addition the counterterm diagram of figure 3, which receives contributions 

from the tr{AiAi) and p'^tr{AiAi) counterterms. These give: 

6 M 16 7 6 4^ 

G,t = -12Zo - 48Zi = -12Zo - 4 H — ) - - 4 + A H^)- (3.65) 

TV jji o TV oTV 

The final result is 



Gcutoff = Gi.ioop+G,t = -12Zo-24M2Ci+8M2C2-^F2+^ \n{^)-^a^) + ^+^. 

(3.66) 



Demanding that this equals the dimensionally regularized result ( 3.60 ) above then 
determines the coefficient Zq to be 

Zo = -— F2 + ln(4^) + -K, - 2M^Ci + -M^C2. (3.67) 
" 15 307r2 187r2 3 ^ ' 

Thus, the required counterterms for our two-loop calculation are given in ( |2.13| ), with 
the coefficients Zq, Zi, and Z2 given in (|3.67| ) and ( |3.18| ), respectively. 
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2a 2b 2c ct 



Figure 4: Diagrams that contribute to the 2-loop free energy. The vertex of 
diagram 2a (2c) corresponds to the first (second) interaction term of ( pj.lOl ). 

4. Evaluating the two-loop diagrams 

Now that we have worked out the required counterterms, we are ready to calculate the 
necessary correlators in (|2.20|) to determine the partition function. The relevant diagrams 



are depicted in figure 4. We begin first with the planar contributions. 

4-1- Two-loop planar diagrams 

From the discussion in section 2.4, the contribution to the partition function from 
planar diagrams comes from the coefficients of tr(t/"^)tr(t/^") in the expansion in powers 
of U. We have three two-loop planar diagrams, one from a pair of cubic vertices and one 
from each of the two quartic vertices. In addition, we have at the same order one-loop 
planar diagrams for each of the counterterms ( |2.13| ). 

We begin with the two-loop diagrams, including a regulator function for each Ai line, 
though we can neglect regulator functions for momenta whose sums are already exponen- 
tially suppressed. Full expressions for these diagrams were previously derived in . They 
arellil (with summation over the spherical harmonic indices a, (3 and 7 implied) 



St = -^^^{D^P'^D^'P^ - I?°^^I?^^^)A,jO,a„,)A,,(0,a„,) 
^g^ia(io + + 2)A,J0, a„6)A,,(0, a„e), 



(4.1) 



^Pj ^ _MM^a/3^^5^7| rftA,Jt,a„,)A,,(t,a,,)A,^(t,a, 



dtAj^{t, aabl^jf, (t, abc)Aj^ (t, a^a), 

(4.2) 



Recall that in diagram 2c, we ignore the part proportional to 5(0). 
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St = P9ym ■ iD^A^^{0,aac)DrA^^iO,aat) + Up + 1)2A, JO, a5c)A,, (0, 



5 Q^6c)Aj^ (0, aab)): 

(4.3) 

where E and the functions i?3± and i?4± are defined in appendix B. Here, the first fine 
for each diagram gives the expression for the diagram before performing any sums over 
angular momenta, while the second line gives the result after summing over everything but 
total angular momentum quantum numbers (these come into the regulator, which we have 
not yet written explicitly). Note that for all diagrams, each of the propagators contributes 
factors of a to two of the three index loops, which we label by a, 6, and c. The notation 
aab indicates that for the tensor products (a ® 1 — 1 (8> a) appearing in the propagator 
( p.l2|) , the first and second elements of the tensor product appear in the traces associated 
with index loops a and 6, respectively. 

We now expand these expressions in powers of U to read off the required coefficients 
fnix) defined in ( |2.21| ), inserting explicit expressions for the regulators at this stage. For 
all diagrams, we find that 

/n(x) = /i(x'^), (4.4) 

so we need only give the results for fi{x). 
We find that diagram 2a gives 

l^a=l a=l ) 

where 

/w-|f^f-'^'-(r^+Mi-.). (4.6) 

For diagram 2b, we find 

/f(^) = f;f; ^"^E^' (6+l)x«+^+2-(a + c+2)x^+i 



2(a + 1)(6 + l)(c + l)(a + 6 + c + 3)(a + c- 6 + 1) 

■ {i2|+(a, 6, c) (a + 6 + c + 3)^ + Rl_ (a, 6; c) (a + 6 - c + 1)^ 

+ Rl_{h, c; a)(6 + c - a + 1)2 + Rl_{c, a; 6)(c + a - 6 + l)^} i?(^)i?(^). 

(4.7) 
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To see the regulator dependent pieces explicitly, we may expand R{^^) about c = a for 
large a, 

_ ,a + l. c-a , a+1 [c - a)^ a + 1 

^^Iw"^ - ^^^^ + -JF^ + ^M^^ + • • • ^^-^^ 

For large a, the sum over c of the summand above with various powers of (c — a) inserted 
then has the following asymptotic behavior as a function of a: 

where we have made the replacements a— s>a— 1, b ^ b — 1. The omitted terms lead only 
to regulator-independent finite contributions. The regulator-dependent terms are then 

fe=l a=l 

_ J_ (b^-l)(86^-7) , 1 .2. A^ 

6=1 a=l 

b=l a=l 

Note that the term involving l/M'Y^^R(a/M)R\a/M) is regulator-independent 
For diagram 2c, we find 

oo oo (a+b)/2 , / L , 1 I 1 \ 

A^- E E E ^^-^^^ - ^ - ^^^^ 



(4.10) 



•{i?2^(a,c,6) + i?t(a,c,6)}i?(^). 



(4.11) 



The asymptotic behavior of the sum over c for large a is 

;2 1 11 /a2 ^\IQU2 



Stt^ b IStt^ a 

The regulator-dependent terms are then 

1 °° 7,2 _ 1 o° ^ 

6=1 a=l 



1 y^ (fe^-l)(86^ + 3) ,^1 g 

6=1 a=l 



(4.13) 
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Combining all the regulator-dependent terms and carefully comparing the divergent sums 
to the divergent integrals defined in (|3.14|) , we find 



= fix)i4C^M'-^C2M'-J^) + ^i8gix)-7fix))H^^) + ^^^^ 

(4.14) 



where f{x) is as above, and 
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(x) = ^b{b'^ - l)x" 



b=i 



(l-x) 



4 ■ 



(4.15) 



The expression ( [4.14]) will be useful in verifying that all regulator dependence cancels, 
though in practice, we simply need to use ( ^75|) , ([4.7|), and (|4.11|) with any convenient 
choice of regulator. 

4-2. Counterterms 

To the expressions above, we must add the counterterm contributions, coming from 



the single-trace counterterms in ( p.l3| ). Using this expression, we find 

Seff = (Set) = f3g^MNj2(^o + (Zi - Z2){jo. + l)')A,jO,a„,). (4.16) 

a 

To determine the term proportional to Z2, we have used the fact that 



-Dl + {j + lf)^,{t,a)=5{t) 



(4.17) 



ignoring the resulting (5(0) term, which will be cancelled by another diagram which is 
independent of U and the temperature. 

Expanding this to find the coefficient of tr(C/"^)tr(C/^'^) and performing the angular 
momentum sum, we find that the counterterm contribution to is /^*(x) = /{^*(x"), 
where 

fl\x) = 2Zof{x) + 2(Zi - Z2)9{x). 



(4.18) 



Using the expressions ( |3.18D and ( |3.67| ) above, this becomes 



MV4C, + ic.) + ^.n(^ 



— F J 

15 ^ 97r2 



+ g{x) 



2^ J 8_ fAi 

15 47r2 157r2 1 ^2 



(4.19) 
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Comparing this with the regulator-dependent pieces ( [4.14|) above, we see that aU regulator 
dependence cancels in the sum " + fi^ + + /f*. Thus, our final result is 

/.(x) = /i(x^), (4.20) 

where is the sum of ( ^4.19|) , (|4.5|) , ( ^4.7|) , and ( [4.11| ), which is finite and regulator 

independent. This may be computed numerically for any value of x, and we do indeed 
find that the numerical results are independent of the regulator used. The perturbative 
expansion of /i at small x takes the form 

14 55 
= —x^ + -x^ + ^x^ + ■ ■ ■ , (4.21) 



and its value at the critical temperature (p..2|) is given by 

/i(xe,o) ^ 0.0253. (4.22) 

The functions fn{x) determine the planar contribution to the two-loop partition function 
via ( |2.23| ) while fi{x) determines the perturbative shift in the Hagedorn temperature, 
which we compute using (|2.26|) in section 4.5 below. 



4-3. Two-loop non-planar diagrams 

We now turn to the non-planar contribution to the two-loop partition function. In 
this case, since there is only a single index loop, and since each term in the propagators 
contributes an equal number of f/'s and C/^'s, we will always end up with just the iden- 
tity matrix inside the single trace. Equivalently, we can simply set a = (C/ = 1) in 
all propagators from the start. The resulting temperature-dependent but t/-independent 
expressions contribute directly to the two-loop free energy via ( p. 24 ). 



We find that the expressions for the three non-planar diagrams are related to the 
expressions (|4.1| )-( |47^ ) for the planar diagrams by setting a = and including an overall 
factor of — Thus, we have 

-2 - _ - 



gnp ^ ^dxM^D^f^JD"^! _ L'°°^L'^^^)AjjO,0)Aj-^(0,0) 

(4.23) 

'^|^Ja(jo + 2)jf,Up + 2)A,„ (0, 0)A,,(0, 0), 



This is related to the fact that the planar and non-planar diagrams must cancel for the 
Abelian theory with A'^ = 1, which is free. 
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dtA.Jt, 0)A,^(t,0)A,^(t, 0), 



(4.24) 



^2T = -/3^7yM— ^-^(^rA,jO,0)I).A,,(0,0) + (j^ + 1)2A,„ (0, 0)A,,(0, 0)) 



+ 2) 

(I^.A,,(0,0)I^.A,„(0,0) + {jp + 1)2A,J0,0)A,,(0,0)). 

(4.25) 

In each case, the group indices are contracted into a single (trivial) trace which gives 
an overall factor of A^. We will find again that the angular momentum sums here are 
divergent, but there are additional counterterms coming from the double-trace terms in 
( p.l3| ). The result for these is obtained from the result for the single-trace counterterms 
in the same way that the non-planar two-loop contributions are obtained from the planar 
two-loop contributions. 

Using this relation between planar diagrams / single-trace counterterms and non- planar 
diagrams/double-trace counterterms, it follows from ( |2.21| ) that the full contribution to 
'^^loop be written as 

oo 

S;f+n, = Pg'vMN J2 /n(x)(tr(C/-)tr(C/tn) _ 1) 

'"=1 (4.26) 
+ P9YM Yl /nm(x){tr(C/^)tr([/-)tr(t/-™) + c.c.-2iV}. 

n>m>0 

For the planar diagrams, the two-trace terms in the first line arose from the combina- 
tion of divergent two-loop contributions plus counterterms to yield the finite regulator- 
independent results defined by fn{x) = /i(x"). While we did not calculate the three trace 
terms explicitly, it is straightforward to show that they are manifestly finite. It immedi- 
ately follows that all divergences and regulator-dependence cancel also for the non-planar 
diagrams plus double-trace counterterms. 

Focusing back on the non-planar terms that we need to calculate, it is convenient to 
use (|4.26|) to write 

oo 

F7{x) = - E + ^''''(^) ' (4.27) 

n=l 
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since we have already computed fi{x); it remains only to calculate 

F7(x) = -2 /-(^)' (4.28) 

n>m>0 

where fnm{x) may be computed from the expressions (|4.1| )-(^7^). Diagram by diagram, 
we find 

1 ia(ia+2)j^(j^ + 2) 



CX3 OO 



(1 - - x^0+^) ~ (1 -a;J/3+i) ~ (1 - ~ ^1 - ^ 

(4.29) 

CX3 OO (ja+j/3-l)/2 

^.T = E E E 

ic = l j/3 = li^/2=(|i„-j0| + l)/2 

{Rl+{ja,jp,jj){ja+jf3+j'y + 3)^ + ^Rl_{ja,jp;j^){jc+jp-j^ + 1)^) 

12(ja + l)(i/3 + l)0'7 + l) 
1 



Ja+jfB+j-y + 3 (1 - XJ« + 1)(1 - XJ'3 + 1)(1 - X-?'^ + l) 

: : + : : + : : + : : - 

Oa + j(3 + 1)X^^ + ^ Ul + l)x^«+^''3+2 



+ — 
+ 



2 / Uf3 + h + l)a;^"+^ (ia + l)x^'5+^'^+2 



(ja + j/3 + + 3) {ip + - + 1) I 1 - XJ" + 1 1 - xJ/3+^^^+2 



(ia + J> + h + 3) (i7 + ia - + 1) I 1 - X^/^ + l 1 - X^'^+J«+2 

(4.30) 

and 

00 00 (j«+J/3)/2 ^ 

^2T=1]1] A {j I 2) J7, i/j) + Rl-UcJ^, jp)) 

jc = lj0 = lj-,/2={\\j^-jp\-l\ + l)/2-^^^'^'^ > 

J/3 + 1 Al+X^«+l)(l + X^'9 + l) (l + X^f + l) (l+X-^° + l) 



+ 1 V(l - - ^^"^^) (l-xJ"+^) (l-xJ°+i) 



+ 1 y (1 - X^«+-?'3+2) 

(4.31) 

Using ( [4.27|) , these results, combined with our results for /i(x) in section 4.2, give the final 
result for -F^^, which in turn gives the non-planar contribution to the two-loop partition 
function and free energy via ( p.23| ) and (|2.24| ). 
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4-4- Results: Two-loop partition function 

We have now calculated all the elements necessary to give the full two-loop partition 
function. Using (|2.23[ ) and ( [4.27| ) , we find that the two-loop partition function for U {N) 
Yang-Mills theory at large on a small is 

~ °° A/3 fi (x") 

where z{x) is given in ( p.l6|) , F2^{x) is the sum of (^4.29| ), ( [4.30| ), and ([4.31j ), and fi{x) is 
the sum of (U), (pi] ), and ( pJ|) . To get the SU{N) result, we simply divide by 

Zu(i)i yielding 

For either U{N) or SU{N), the order A correction to the free energy (using (|2.24| )) is 



SF = -Td ln{Z) = A 



F-w + (^3^-1) 



(4.34) 



n=l 

Our results for the two-loop partition function can be expanded in powers of x as in ( |1 . 1|) , 
from which we can read off the sum of energy corrections for all states with a given energy 
in the free theory, since 

J2 x^o+x5E^ ^ ^Eo ln(x) + ^(-^') j • (4-35) 

For SU{N), we find that the first few terms in the series give 

zSo^'ol = 1 + (^21 + ^ A ln(x)^ x^ + (^96 + ^ A ln(x)^ + (^392 + A ln(x)^ x^ + ---. 

(4.36) 

Thus, for example, the sum of perturbative energy shifts for the 21 independent states 
with energy 4/i?53 is 4A/7r^. 

Since the free Yang-Mills theory in the small volume limit is conformal, we have a 
map between states in this theory and local operators in Euclidean Yang-Mills theory on 
IR^. The coefficient of the term proportional to x"^ln(x) is then interpreted as the sum of 
anomalous dimensions for all operators with dimension n. Thus, we can use the results 
of [|2| for anomalous dimensions of operators in pure Yang-Mills theory as a check on our 
results (see section 5.1 below). 
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4-5. Results: Order A correction to critical temperature 

As explained in section 2.4, the two-loop corrections to the partition function will shift 
the critical temperature associated with the Hagedorn and deconfinement transitions. To 
compute the shift in Xc or T^,, we use (|2.26|) or (|2.27|) , which only involve the result ( ^4.22| ) 
for /i(xc). Recalling from section 2.3 that Xc,o = 2 — y/S for A = 0, we find 

5xc = 0.00298A, (4.37) 

so the transition occurs at 

Xc = 2- v^+ 0.00298A + C(A2). (4.38) 
Equivalently, from ( p.27|) , the critical temperature to order A is 

^ r.,„«. . (l . . 0(A^)) . (4.39) 

Below, we will reproduce this result using a formula in ^ for the shift in the transition 
temperature from anomalous dimensions, providing a strong quantitative check of our 
results. 



5. Checks 

5.1. Anomalous dimensions for low-dimension operators 

As a first check on our results, we will try to reproduce our result ( [4.36| ) for the 
leading terms in the expansion of the SU{N) partition function using the results of for 
anomalous dimensions of operators in pure Yang-Mills theory. 

Explicit results for the one-loop anomalous dimensions for low-dimension operators 
in pure large SU{N) Yang-Mills theory are given in table 1 of 0. They show that 
two of the dimension four states have anomalous dimension — (ll/3)(A/87r^), ten have 
anomalous dimension (7/3)(A/87r^), and nine have vanishing anomalous dimension, while 
at dimension five the 16 primary states have anomalous dimension 3(A/87r^) (and each of 
the dimension four primary states that has an anomalous dimension has four descendants) . 
Adding also all the dimension six states, we find that the results of lead to 

Z2-ioop{x) = 1 + (^21 + ^ ln(x)^ + + ^ ln(x)^ x^ + (^392 + ^ ln(x)^ + ■ ■ ■ . 

(5.1) 

This precisely agrees with our result (|4.36|) , if we recall that, due to a difference in con- 
ventions, A in is twice the A that we used in our computation (as was also the case for 
the comparison of the flat-space two-loop partition functions in [H). 
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5.2. Shift in Tc from anomalous dimensions 

In this section, we use an alternate method, based on [Q, to check our result for the 
order A shift in the transition temperature. 

Single trace operators in pure Yang-Mills theory may be put in one to one corre- 
spondence with spin chains of lengths / = 1,2, ■■■,00, where the spins in these chains are 
vectors in one of two representations of the conformal group. The primaries of these rep- 
resentations are the self-dual and anti-self-dual parts of F^j^ and carry quantum numbers 
(1, 0, 2) and (0, 1, 2), respectively, under (ji, J27 A) {ji labels the representation of the first 
SU{2), j2 of the second SU{2), and A is the scaling dimension). We will call these two 
representations chiral c and anti-chiral c, respectively. It was demonstrated in that 

00 

c X c = (0, 0, 4) + (1, 0, 4) + (2, 0, 4) + ^(2 + 1,1,4 + j). (5.2) 
The tensor product of c with c is obtained from ( |5.2| ) by a ji ^ j2 flip. Finally, 



cxc = 5^(|,|,2-Fj)- (5.3) 

J =2 

We now list xr = T^Ri^^)^ the characters for all these representations (they will be 
needed below). Actually, we will flnd it more convenient to list x'r = (1 ~ x)'^XR '■ 

Xo,0,4 ~ ' -^0,1,4 ~ "^-^ ' Xo,2,4 ~ '^•^ ' 

X2+i,|,4+j = (5+i)(i + l)x4+^' -i(4+i)x5+^ (5.4) 

According to equation (6.10) of 

bx^ _ Aln(xc,o) , , , . 

where z{x) was deflned in ( p.l6|) , {D2{x)) is (see equations (4.4) and (4.8) in [|1)0 
11 1 7 00 , 

{D2{x)) = - yX0,0,4 + 2X0,1,4 + 2X0,2,4 + Yl [^^(^ ^ ~ y ) ^2+i,i,4+j 

(5-6) 

+ ^(hU - 2) + h{j + 2) - ll/6)x,,,,2+„ 
J =2 



Note that here, unlike the previous subsection, we are using the same conventions for A as 
in the rest of the paper. 
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and h{j) = ^{^^ 1/k (with h{0) = 0). 

It is easy to numerically evaluate the sum (|5.6| ) at the Hagedorn temperature, or to 
compute it analytically using 



{D2ix)) 



1 



(1-X)6 

Using (|5.5[) , we find 



2(- ln(l - x))(l - 3x - + x^f - 2x + llx^ - -x^ - 26x* + 36x^ 

(5.7) 



{D2{xc,o)) = (28 - leVs) / (V3 - l)' = 1 ^ STjT.^o = j^{D2{x,,o)) 



A 



(5.8) 



or 



T,i?s3 = T,,oi?s3 ■ (1 + + 0(A 



127r-^ 



Finally, given 



(2 - In (2 + VS) 



(2 - 73) In (2 + V3) 
12^ 



(2 - VS) In (2 + v^) a' 



127r2 



0.00298 



(5.9) 



(5.10) 



we find perfect agreement with ( [4.39| ) (within our numerical accuracy). Note that, in 
general, in order to have agreement between (|2.26| ) and (^]^), we should have 

{D2{Xc,0)) 



fl{Xc,o) 



47r^ 



(5.11) 



which is indeed obeyed (within our accuracy) by our result (|4.22| ). 
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Appendix A. Actions and propagators 



The quadratic action for gauge-fixed Euclidean pure Yang-Mills theory on x S^, 
as described in §2, is given by (with Dj^ = dt — i[a, ■]) 

S2 = J tr(^A"(-D^ ^ ^ ^^2^^a ^ ^a^'j^U^ + 2)a" + c^jaOa + 2)c"). (A.l) 
In addition, we have cubic interactions 



= 9YM J c/t tr(zc"[A^, c^]C"^'^ + 2za"A^a^C"^^ 

and quartic interactions 

^4 = 9ym J dttv{ - A%\ A'] (^D^'^^D'^^ -f -^-L-^C«^^C^^^^ 

_ -A°'A^A^A^ ^]J'^l>']jPS\ _ jjaSXjjlS-yX^y 



(A.2) 



(A.3) 



The quantities C, D, and E are integrals over spherical harmonics, and are defined in 
appendix B. The propagators of the various fields follow from ( [A.l| ) and are given by 

ic^Bity^t')) = ■ ,Ao/ "^^(^ - t')Sad5cb. (A.4) 

(«:.(^)«f.(^')) = -TTTV.S-^^ii - i')SadSct. (A.5) 

(^a.(t)^f.(t')) = 5"^A;f '^'(t - t\ a), (A.6) 
where A is defined in section 2. 



Appendix B. Spherical harmonics on 

A detailed discussion of spherical harmonics on may be found in appendix B of |P] . 
Here, we collect various formulae relevant for the present calculation. Many of the basic 



results were derived in [12 . 
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B.l. Basic properties of spherical harmonics 

Scalar functions on the sphere may be expanded in a complete set of spherical har- 
monics SJ^"^ transforming in the (j/2,j/2) representation of SU{2) x SU{2) = SO(4), 
where j is any non-negative integer, and < m,m' < j/2. It is convenient to denote 

the full set of indices (j, m, m') by a. These obey an orthonormality condition (we take 
the radius of the to be one) 



I S3 

where 5" denotes the complex conjugate of S"", 

The spherical harmonics are eigenfunctions of the Laplace operator on the sphere, 

V25" = -ia(ja + 2)5", (B.3) 

and under a parity operation transform with eigenvalue (—I)-'". 

A general vector field on the sphere may be expanded as a combination of gradients 
of the scalar spherical harmonics plus a set of vector spherical harmonics VJI^. . These 
transform in the (^^^, '^-j^) representation of S'0(4), where j is a positive integer. Again, 
it is convenient to denote the full set of indices (j, m, m', e) by a single index a. These 
obey orthonormality relations 

(B.4) 

y° . = . 

Again indicates the complex conjugate of V", given by 

The vector spherical harmonics are eigenfunctions of parity with eigenvalue (—1)-'"'"^, and 
satisfy 

Vxy- = -e,(j,-M)y", (B.6) 
V ■ = 0. 

Explicit expressions for the scalar and vector spherical harmonics may be found in 

li- 



as 



B.2. Spherical harmonic integrals 

The vertices of the mode-expanded Yang-Mills theory on have coefficients involving 
integrals over three spherical harmonics. We define 

JS3 



S3 



S3 



It is also convenient to define 



iaOa + 2)i/3(i/3 + 2) J 33 
^ 1 OaOa + 2) + jf3{j(3 + 2) - j^jj^ + 2)) ^^^^ 

2 i„(ia + 2)i^(i^ + 2) 

^0/37 ^ E^P\t^{3^ + 1) + e;3(i^ + 1) + e^(i^ + 1)) . 



These integrals were calculated in [jT2|, and the results may be expressed asll 



la J_§_ 3j_ 

2 2 2 



where 



ttIq mp rrij / \rn^ rrip 

2 2 2 



2£ 2t_ 
2, 2^ 2^ ) RlUo.JpJ'r)^ 



2 2 2 

2 2 2 

rrict mp ruj 







2 




m'p 




Jet ^CK 


i/3-e/3 


2 


2 







3a—<^a 313— 3i—^-i 
2 2 2 

m!^ 



a 



(B.7) 



:B.8) 



(B.9) 



(-1)- /(x-M)(2/-M)(;2-M) 



TT 



R2{x,y,z) 



TT 



(x + + 1)((T^ - x){a' - y){a' - z){a' + 1) 

+ 



(B.IO) 



fB.lll 



The expression for C below differs by a factor of two from the expression in [ |12| |, but we 
believe that this expression is correct. 
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■ ((e^(x + l) + €,{z + l)+y + 2)(e^(x + 1) + €,{z + 1) + y) 
{e^{x + 1) + e^{z + 1) - y){e^{x + 1) + e^{z + 1) - y - 2))\ 

^ (-1)"'+' . , f(a' + l)(a' -x)(a' -y)(a' - z)Y 

4e.6„6,v .y. ; ^ B\x y zjy 4(x + l)(y + l)(z + l) / 

■ ((e,(x + 1) + ej,(y + 1) + e,(z + 1) + 2)(e,(x + 1) + (y + 1) + e,(z + 1) - 2)) 

(B.12) 

Here, the right-hand sides of the equations are defined to be non-zero only if the triangle 
inequality \x — z\<y<x + z holds, and if a = {x + y + z) /2 (in Ri and R^) and 
a' = {x + y + z + 1) /2 (in R2 and R4) are integers. We also define i?3+ = -R3++ = -R3 , 

i?3_ = i?3^ — = i?3 — 1_, i?4+ = i24+_|_+ and i24_ = i?4+_|__. 

iJ.^. Identities for sums of spherical harmonics 

Using identities for 3j-symbols which may be found in 0, it is straightforward to 
derive expressions for sums over m, m', and e in various products of the spherical harmonic 
integrals. For our calculations, we require: 

m's 

V I^"^^ = — 5A,0ja(ia+2), 

m's,e 

E = 2i?l+(j«, j^, j^) + 2Rl_{j^J^Jf,), 

m' s,e' s 

J2 D^P-^D^'P^ = ^UJc + 2)jp{jp + 2), 

m's 

Appendix C. Useful formulae for dimensional regularization 

The following formulae are useful for our calculations in dimensional regularization of 
section 3.2: 



(B.13) 



Cf^k A;2"^ _ 1 1 r(?i - m - f ) r(m + f ) 

(2^(PTm2)^ - M^n-2m-d i^^^^i Y{^^ r(|) ' 



(C.l) 
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C(l + e) 



1 



+ 7 + 



(C.2) 



e 



r'(^) 



(C.3) 




n— 1 ^ 

A;=l 



(C.4) 



^(_+n) = _^_21n(2) + 25^ 



1 



(C.5) 



2k -1 



k=i 



Appendix D. Regulating the Coulomb gauge 



For our calculations here and in |]5[, we have used the Coulomb gauge, setting the 



divergence of the spatial gauge field to zero. In this gauge, the time component of the 
gauge field and the ghosts have no kinetic term, so their propagator contains a delta 
function. For diagrams containing Aq or ghost loops, we then get a 5(0) factor. While 
this always cancels between the diagrams, one may be concerned that certain finite terms 
have been missed. Further, the singular propagators lead to ambiguities (which we did 
not encounter here) in certain calculations when the time derivative of the Ai propagator 
must be evaluated at t = 0. 

To give an ambiguity-free definition of the Coulomb gauge, we can replace the strict 
V ■ A = condition by a gauge-fixing action 



The Coulomb gauge may be defined as the ^ — > cxd limit of this, but performing calculations 
at finite ^ avoids any 5(0) singularities (they show up as terms that cancel between the 
diagrams). In practice, we need only keep the and terms for each diagram. Using 
this procedure, we can verify that no additional finite terms arise in the cancellation of the 
5(0) 's in our calculation. 

On the other hand, we have found that for certain other calculations, the naive 
Coulomb-gauge calculations can miss finite contributions. As an example, in the two-loop 
Casimir energy on the sphere, a diagram with two jcAqc vertices gives a ^-independent 
contribution that would be missed if we set ^ = oo from the start. 





together with the corresponding ghost action 




(D.2) 
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